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By using the most general form of Einstein equations for General Relativistic (GTR) spherical collapse of an 

. -j ' isolated fluid having arbitrary equation of state and radiation transport properties, we show that they obey a Global 

Constraint, 2GM(r,t)/R(r,t)c 2 < 1, where R is the "invariant circumference radius", t is the comoving time, and 

m : M(r ' t] is the ^ vitational mass enclosed within a comovi ^ sheU '■ This inequality specifkally shows that < contrary 

to the traditional intuitive Newtonian idea, which equates the total gravitational mass (Mb) with the fixed baryonic 

^. . mass (Mo), the trapped surfaces are not allowed in general theory of relativity (GTR), and therefore, for continued 

collapse, the final gravitational mass Mf — > as R — v 0. This result should be valid for all spherical collapse scenarios 

including that of collapse of a spherical homogeneous dust as enunciated by Oppenheimer and Snyder (OS). Since 
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the argument of a logarithmic function cannot be negative, the Eq. (36) of the O-S paper (T ~ In v l _ l ) categorically 
demands that yb = Rb/R g b > 1, or 2GMb/RbC 2 < 1, where Rb referes to the invariant radius at the outer boundary. 
Unfortunately, OS worked with an approximate form of Eq. (36) [Eq. 37], where this fundamental constraint got 
obfuscated. And although OS noted that for a finite value of M(r, t) the spatial metric coefficient for an internal point 
fails to blow up even when the collapse is complete, e A ( r<ri> ) ^ oo for R — ► 0, they, nevertheless, ignored it, 
cZ) ■ and, failed to realize that such a problem was occurring because they were assuming a finite value of Mf, where Mf 

a ■ 

L* ' is the value of the finite gravitational mass, in violation of their Eq. (36). 

Additionally, irrespective of the gravitational collapse problem, by analyzing the properties of the Kruskal transfor- 



mations we show that in order that the actual radial geodesies remain timelike, finite mass Schwarzschild Black Holes 
cannot exist at all. 

Our work shows that as one attempts to arrive at the singularity, R — > 0, the proper radial length I = 
J \J—9rr dr — > oo (even though r and R are finite) , and the collapse process continues indefinitely. During this 
indefinite journey, naturally, the system radiates out all available energy, Q —> MiC 2 , because trapped surfaces are 
not formed. And this categorically shows that GTR is not only "the most beautiful physical theory" , but also, is the 
only, naturally, singularity free theory (atleast for isolated bodies), as intended by its founder, Einstein. However, this 
derivation need not rule out the initial singularity of "big bang" cosmology because the universe may not be treated 
as an "isolated body". 

There is a widespread misconception, that recent astrophysical observations have proved the existence of Black 
Holes. Actually, observations suggest existence of compact objects having masses greater than the upper limit of 
static Neutron Stars. The present work also allows to have such massive compact objects. It is also argued that 



there is evidence that part of the mass-energy accreting onto several stellar mass (binary) compact objects or massive 
Active Galactic Nuclei is getting "lost" , indicating the presence of an Event Horizon. Since, we are showing here that 
the collapse process continues indefinitely with local 3-speed V — > c, accretion onto such Eternally Collapsing Objects 
(ECO) may generate little collisional energy out put. But, in the frame work of existence of static central compact 
objects, this small output of accretion energy would be misinterpreted as an "evidence" for Event Horizons. Thus the 
supposed BHs are actually massive compact ECOs. 
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I. INTRODUCTION 

One of the oldest and most fundamental problems of physics and astrophysics is that of gravitational collapse, 
and, specifically, that of the ultimate fate of a sufficiently massive collapsing body [jjlp). Most of the astrophysical 
objects that we know of, viz. galaxies, stars, White Dwarfs (WD), Neutron Stars (NS), in a broad sense, result 
from gravitational collapse. And in the context of classical General Theory of Relativity (GTR), it is believed that 
the ultimate fate of sufficiently massive bodies is collapse to a Black Hole (BH) ||. A spherical chargeless BH of 
(gravitational) mass Mj is supposed to occupy a region of spacctime which is separated by a hypothetical one-way 
membrane of "radius" R g b — 2GMj ) /c 2 , where G is the Newtonian gravitational constant and c is the speed of light. 
This membrane, called, an event horizon, where (local) acceleration due to gravity blows up, is supposed to contain 
a central singularity at R = 0, where the other physically relevant quantities like (local) energy density, (local) tidal 
acceleration, and components of the Rimmenian curvature tensor diverge. However, although such ideas are, now, 
commonly believed to be elements of ultimate truth, the fact remains that, so far, it has not been possible to obtain 
any analytical solution of GTR collapse equations for a physical fluid endowed with pressure (p), temperature (T) 
and an equation of state (EOS). And the only situation when these equations have been solved (almost) exactly, is by 
setting p = 0, and further by neglecting any density gradient, i.e., by considering p = constant B. It is believed that 
these (exact) asymptotic solutions actually showed the formation of BH in a finite comoving proper time T g b- However, 
this, assumption of perfect homogeneity is a very special case, and, the speed of sound c s = (dp /dp) 1 ' 2 = oo in such 
a case if p > 0. And this is not allowed by GTR or Special Theory of Relativity (STR). In fact, now many authors 
believe that for a more realistic inhomogeneous dust, the results of collapse may be qualitatively different [pH7|. These 
authors, on the strength of their semi- analytical and numerical computations, claim that the resultant singularity 
could be a "naked" one i.e., one for which there is no "event horizon" or atleast some light rays can escape the 
singularity. There has also been some effort to study the final stages of collapse by assuming the presence of positive 
and likely negative pressure gradients subject to the inherent difficultities and limitations of any "direct approach" . 
For instance, recently Cooperstock et al. M undertook such a study and tentatively concluded that for positive 
pessure gradients a BH is likely whereas for occurrence of negative pressure gradients a "naked singularity" may form. 
Therefore light may emanate from a naked singularity and reach a distant observer. A naked singularity may also 
spew out matter apart from light much like the White Holes. In other words, unlike BHs, the naked singularities are 
visible to a distant observer and, if they exist, are of potential astrophysical importance. However, according to a 
celebrated postulate by Penrose M, called "Cosmic Censorship Conjecture", for all realistic gravitational collapse, the 
resultant singularity must be covered by an event horizon, i.e, it must be a BH. And many authors believe that the 
instances of occurrences of "naked singularities" are due to fine tuned artificial choice of initial conditions or because 
of inappropriate handling and interpretation of the semi-analytical treatments. In this paper, we are not interested 
in such issues and would avoid presenting and details about the variants of naked singularities (strong, weak, local, 
global, etc.) or the variants of the censorship conjecture. 



A. Kelvin - Helmholtz (KH) Process 

As a self-gravitating body contracts, it radiates and the same time its internal energy increases. And the internal 
energy can have two contributions: 

E in = Et + E co id (1.1) 

where Et — J epdV is the temperature dependent thermal part of the internal energy and E co id — J e co iddV is due 
to the pure degeneracy effects and which may exist in certain cases even if the star is assumed to be at a temperature 
T = 0. Here dV is an element of proper volume. The corresponding energy densities are 

(3/7r 2 ) 1/3 mc 2 .,,- 

Here m is nucleon mass and n is nucleon number density. Actually, when the body is really degenerate, this kind of 
splitting of E in can be done only in an approximate manner. For example, if it is assumed that a degenerate ideal 
neutron gas is close to T = 0, i.e, T <C Tp erm i, then one may approximately take the first term (lowest order in T) 
of an infinite series to write the above expression. On the other hand, if the temperature is indeed much higher than 
the corresponding Fermi temperature, degeneracy will vanish, e co id — ► 0, and the entire energy density will be given 
by the thermal contribution: 
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e = er = —nkT; T — > oo (1-3) 

where k is the Boltzman constant, and 

2 2(3tt 2 ) 2 / 3 h 2 5/3 

e co id = rPcoid = -^—ri n 5/d ; j t = 5/3 (1.4) 
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where "ft is the effective ratio of specific heats. Since it is not known beforehand, how T would evolve, in principle, 
one should work with an expression for ex (an infinite series) valid for arbitrary T. But, it is not possible to do so 
even for an ideal Fermi gas. As to the actual EOS of nuclear matter at a finite T, it may be remembered that, it is 
an active field of research and still at its infancy. Thus, in practice it is impossible to make much headway without 
making a number of simplifying assumptions because of our inability to self consistently handle: (i) the equation of 
state (EOS) of matter at arbitrarily high density and temperature, (ii) the opacity of nuclear matter at such likely 
unknown extreme conditions, (iii) the associated radiation transfer problem and all other highly nonlinear and coupled 
partial differential (GTR) equations (see later). 

One may start the numerical computation by presuming that indeed the energy liberated in the process Q <C M^c 2 , 
i.e., the effect of GTR is at best modest. Then, it would naturally be found that the temperature rise is moderate and 
depending on the finite grid sizes used in the analysis and limitation of the computing machine, one may conclude 



that the formalism adopted is really satisfying, and then find that Q <C MiC 2 [pLO 11|. Meanwhile, one has to extend 
the presently known (cold) nuclear EOS at much higher densities and maintain the assumption that the rise in 
temperature is moderate. Because if T is indeed high, in the diffusion limit, the emitted energy Q ~ T 4 would be 
very high, and the value of Mf could drop to an alarmingly low value. Thus, for the external spacetime, one needs to 
consider the Vaidya metric ]l^-fi4|. Actually, even when, T is low, it is extremely difficult to self consistently handle 
the coupled energy transport problem. 



It may appear that, the practical difficulties associated with the study of collapse involving densities much higher 
than the nuclear density can be avoided if one starts with a very high value of Mi, say, 1O 1O M (solar mass). Then 
if one retains the assumption that Mi = Mf, one would conclude that an "event horizon" is formed at a density of 
~ 1CP 4 g cm~ 3 for which the EOS of matter is perfectly well known: 

o» = 3d^r 2 * 10,a e cm - 3 (£)" (1 ' 5) 

Here p g is the closure density. It may be reminded that this p g or the generic mass-energy density in GTR p does not 
include the (negative) contribution of self-gravity. On the other hand, p includes rest mass and all internal energy 
densities. It is well known that, in GTR, there is no locally defined gravitational energy density; however the effect 
of (negative) self-gravity appears in globally defined concepts like enclosed Schwarzschild (gravitational) mass: 
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M = / pdV (1.6) 

Here dV is an element of "coordinate volume" element and not the locally measured physically meaningful "proper 
volume" element dV. While p refers only to the material energy density and does not take into account the contribution 
to the energy from gravity and while dV is the coordinate volume and not the physical proper volume, dV ', the 
combination of these elements in Eq. (1.6) work together to yield the correct total energy of the body (including) the 
contribution from gravity [15] . What is overlooked in the traditional interpretation of Eq. (1.5) is that, this expression 
is incorrect, and the correct expression should involve Mf and not M^. 
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Once we are assuming that an event horizon is about to form, we are endorsing the fact that we are in the regime of 
extremely strong gravity, and, therefore for all the quantities involved in the problem, a real GTR estimate has to be 
made without making any prior Newtonian approximation. 

To further appreciate this important but conveniently overlooked point (by the numerical relativists), note that, 
the strength of the gravity may be approximately indicated by the "surface redshift" , z s , of the collapsing object, and 
while a Supermassive Star may have an initial value of z s as small as 10 -10 , a canonical NS has z s ~ 0.1, while the 
Event Horizon, irrespective of the initial conditions of the collapse, has got z s — oo! Therefore all Newtonian or Post 
Newtonian estimates or the conclusions based on such estimates have little relevance for actual gravitational collapse 
problem. 

As a result, the integrated value of Q may tend to increase drastically, and this would pull down the running value 
of Mf = Mq — Q/c 2 and R g b to an alarming level! At the same time, of course, the value of R is decreasing. But 
how would the value of Rb/R g b would evolve in this limit? Unfortunately, nobody has ever, atleast in the published 
literature, tried to look at the problem in the way it has been unfolded above. On the other hand, in Newtonian 
notion, the value of Mf is permanently pegged at Mq because energy has no mass-equivalence (although in the 
corpuscular theory of light this is not so, but then nobody dragged the physics to the R — > R g limit seriously then). 
So, in Newtonian physics [jj]g], or in the intuitive thinking process of even the GTR experts [|J, the value of 

2GM 2GM , s 



and the idea of a trapped surface seems to be most natural. But, in GTR, we cannot say so with absolute confidence 
even if we start with an arbitrary high value of Mj because, in the immediate vicinity of R — > R gi the running value 
of M may decrease in a fashion which we are not able to fathom either by our crude qualitative arguments, based 
on GTR, or by numerical computations plagued with uncertain physics and inevitable machine limitations. And, if 
Mf drops to an alarming level, the actual value of p g can rise to very high values. Thus all the difficulties associated 
with the numerical study of the collapse of a stellar mass object may reappear for any value of Mq unless one hides 
the nuances of GTR and other detailed physics with favorable and simplifying assumptions and approximations. To 
seek a real answer for such questions, we need to handle GTR carefully and exactly in a manner different from this 
qualitative approach. 

II. SPHERICALLY SYMMETRIC GRAVITATIONAL FIELD 

The most general form of a spherically symmetric metric, after appropriate coordinate transformations, can be 
brought to a specific Gaussian form [I15H18| , : 

ds 2 = A 2 (xi,x )dxl - B 2 (x Q ,x 1 )dx 2 1 - R 2 (d9 2 + sin 2 9d<\> 2 ) (2.1) 

where A and B are to be determined self consistently for a given problem. For the simplification of computation, it 
is customary to express: 

A 2 = e u - B 2 = e A (2.2) 

Here x\ is an appropriate radial marker (coordinate), and xq is the coordinate time. For the spherical metric, if we 
consider aio = constant hypersurface and pick up a curve (circle) with x± — constant and 9 = n/2, the value of the 
invariant line element would be 

ds = Rdcf) (2.3) 

The invariant circumference of the xq = constant circle would be 2ttR, and thus, we identify R as the (invariant) 
circumference variable. The R thus defined acquires a physical significance and appears to be related to the luminosity 
distance for astronomical observations. Clearly this property of R as the invariant circumference radius is a pure 
consequence of spherical symmetry irrespective of whether there is a probable coordinate singularity or not. One has 
a Schwarzschild coordinate system when one chooses x\ = R. If we label the corresponding time as T, we have 

ds 2 = e u dT 2 - e x dR 2 - R 2 (d6 2 + sin 2 9d<j) 2 ) (2.4) 

And if we demand that e v — ► 1 as R —> oo, we can identify T as the proper time of a distant observer. 

A. Exterior Schwarzschild Metric (ESM) 

But, now, suppose we are going to describe a truly "vacuum" exterior spherical solution, an exterior spacetime 
region not containing a single "particle" or photon. The actual solution for the vacuum exterior spacetime region was 
found by Schwarzschild in 1916 [ |19[ : 
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ds 2 = g T rdT 2 + g RR dR 2 - R 2 (d9 2 + sin 2 9d(j) 2 ) (2.5) 

with 

9tt = (l- -|-J ; g RR = -[i- -|-J ; R gb = —^- (2.6) 

where Mb is the total gravitational mass of the system. The time parameter T appearing here is not any comoving 
time measured by a clock attached to the test particle. On the other hand, as mentioned above, T acquires a distinct 
physical meaning as the proper time measured by a distant inertial observer Soo- Thus both R and T have some sort 
of absolute meaning in the External Schwarzschild Metric (ESM) . The singularity of this vacuum mertic at R = R g b 
is obvious, and, as per Landau and Lifshitz |]15| , if any fluid is squeezed to R < R g b, this singularity means that, the 
body cannot remain static in such a case. On the other hand, the body must be collapsing, if it ever reaches the 
Schwarzschild surface. If during the preceding collapse process, the body radiates, the vacuum condition would break 
down and the metric would simply become inappropriate to describe a radiating scenario. However, if the fluid is a 
dust, the collapse process would be radiation free, and the vacuum solution must continue to hold good exterior to 
the collapsing dust at every stage because the derivation of Eq. (2.6) is absolutely general (except of course for the 
"vacuum" condition) . 

In contrast, a comoving frame (COF), by definition, can be constructed in a region filled with mass-energy and 
can be naturally defined in the interior of any fluid. And since in this case x\ is fixed with the particle and time is 
measured by the clock moving with the test particle, the question of any coordinate singularity does not arise. 

To appreciate, again, the point that R always retains it physical significance as the "metric distance" for spherical 
symmetry, consider the case of the infall of a free particle towards the central singularity of a supposed massive BH 
(if it would exist). The conventional wisdom is that when the particle crosses the Event Horizon, R and T would 
exchange their roles, i.e, R would be time like while T would be space like. However, note that, even in this case the 
"location" of the central singularity is still denoted by R = and not by T — or any T = Tq. Further, if one would 
like to introduce a totally different coordinate system, such as Kruskal Coordinates u and v (see later), the central 
singularity is still expressed, most conveniently, as R = and not as u = or v = 0. Even if one would mechanically 
express the central singularity as, say, u = uo, even then, the idea of R — stalks in the back ground, and u needs to 
be calibrated against R for a meaningful physical description. Without the support of the concept of R, the essential 
spatial distinction between two events cannot at all be described. 

III. FORMULATION OF THE COLLAPSE PROBLEM 

The general formulation for the GTR collapse of a perfect fluid, by ignoring any emission of radiation, i.e, for 
adiabatic collapse, is well developed ]lq , p0| , pi[ . Although, our central result would not depend on the details of the 
numerous equations involved in the GTR collapse problem, yet, for the sake of better appreciation by the reader, we 
shall outline the general formulation of the GTR spherical collapse problem, and refer, the reader to the respective 
original papers for greater detail. It is most appropriate and simple to formulate the collapse problem in the comoving 
frame xq — t and x\ = r: 

ds 2 = A 2 (r,t)dt 2 - B 2 (r,t)dr 2 - R 2 (r, t)(d6 2 + sin 2 6dcj) 2 ) (3.1) 



The most physically significant choice for r is one which corresponds to the fixed number of baryons inside a shell of 
r = r. The element of proper time is naturally obtained by following the particle r = constant [ p.5| in radial motion 

ds = dr = Adt (3.2) 

Now we set c = 1. On the other hand, for this static metric, the element of proper length along the radial motion of 
the fluid is |0§ 

dl = Bdr (3.3) 

At first sight, one might think if r is really a "comoving coordinate" and dr — for a particular r, then what does 
a derivative with respect to r or a derivative for r mean? Here dr is to be interpreted as the difference in the value 
of r between two close by shells having r — r and r = r + dr, and also the increment of a particular shell during the 
overall fluid motion. To understand the latter interpretation of dr, we can first consider a Newtonian hydrodynamics 
problem. Suppose we take snapshots of the fluid shells corresponding to comoving times t and t + dt recorded by a 
given clock fixed initially at r = r. Now, if we superimpose the two snapshots, we can measure the the radial increment 
of a particular shell by comparing its relative position against the backdrop of the other. The corresponding derivative 
dr/dt will give the local speed of the fluid in terms of comoving coordinates in the Newtonian case. In the GTR case 
too, one can similarly define a derivative dr/dt and use it for further defining appropriate 3 -speed (see later). Before 
we proceed, we may define few variables which, to start with may be taken as pure symbols [fl8|.p0| |23|]: 



and 



tt dR dR R 

U ^^ = Adt l- = A (3 - 4) 



dR dR R' 

T = ^f=Bd-r lt = t= B (3 - 5) 

where a prime denotes differentiation with respect to r and an overdot denotes the same with respect to t. Note that 
while U has the nature of a partial derivative with respect to t, it is a total derivative with respect to r because the 
concept of a fixed r is inherent in the concept of dr. Similarly, while T has the nature of a partial derivative with 
respect to r, it is a total derivative with respect to I because the concept of a fixed t is inherent in the concept of 
element of proper length dl (see pp. 180-181 of ref. 22 and pp. 150-151 of ref. 23). 
We now recall the Einstein equation itself: 

Rik = 8ttG (r ik - ~t) (3.6) 

where the energy momentum stress tensor for a perfect fluid, in the COF, is 

T r r =T e e =T* = ~p; T °=p; T = T\ = p - 3p (3.7) 

Here p is the isotropic pressure (in the proper frame) and the total energy density of the fluid in the same frame 
(excluding any contribution from global self-gravitational energy) is 

p = Pa + e (3.8) 



where po = mr >' is the proper density of the rest mass, n is the number density of the baryons in the same frame (one 
can add leptons too) , and e is the proper internal energy density. Here Rik is the contracted (fourth rank) Rimennian 
curvature tensor R^ M , and, is called the Ricci tensor. In terms of the Christoffel symbols 
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the components of the Ricci tensor are 
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dx l dx k 
One also requires to use the local energy momentum conservation law 

where a semicolon, ";" , denotes covariant differentiation: 
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One has to supplement these equations with the equation for continuity of baryon number : 

(nu l ):i = 
Now, after considerable algebra, in the COF, the Einstein equations become [El| : 

«): ^G P R 2 R>= 1 -(r+^-^-) 



(3.9) 
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(3.15) 



(R 6 e ,Rt): A7rG(p + p)R 3 =[R 
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(3.16) 



(Rq, R r ) : - — — - + — R' 

Further, if we define a new function 

M(r, t) = 4?r / pR 2 R'dr 
Jo 

the Rq- field Eq. (3.14) can be readily integrated to 

RR 2 RR' 2 



R 



= 2GM 



(3.17) 



(3.18) 



(3.19) 



A 2 B 2 

where the constant of integration has been set to zero because of the standard central boundary condition M (0, t) = 0. 
Here the coordinate volume element of the fluid is dV = AnR-R'dr. If we move to the outermost boundary of the fluid 
situated at a fixed r — rj, and demand that the resultant solutions match with the exterior Schwarzschild solution, 



then we would be able to identify M(rb) as the function describing the total (gravitational) mass of the fluid as 
measured by the distant inertial observer Soo. For an interpretation of M(r,t) for the interior regions, we first, recall 
that the element of proper volume is 

dV = 4nR 2 dl = 4irR 2 Bdr = — (3.20) 

so that 

pT pr 

M(r,t)= / pdV = / [Tp)dV (3.21) 

Jo Jo 

And this suggests that Tp is the energy density measured by S^. We have already inferred that Mf, is the total 
mass energy of the fluid as seen by Soo. Then for a self-consistent overall description, we can interpret that, in general, 
M(r, t) is the mass-energy within r = r and as sensed by Soo fl20|| . 

Now by using the definitions U and T from Eqs. (3.4) and (3.5) in Eq. (3.19), we find that 

r 2 = 1 + U 2 — (3.22) 

R 



Id, Id, , 

Dt = Adt^ ^ = --\ t , t (3.23) 



Further using a compact notation 



the major adiabatic collapse equations turn out to be [ p0|j21[ | : 

D T M = 4irR 2 pT (3.24) 

D t M = -4irR 2 pU (3.25) 

ftp— ^Lf*! (3.27) 



p + p \dR 



An immediate consequence of the last equation is that, if we assume a p = EOS, T will be time independent 
r(r, t) = r(r), and for a fixed comoving coordinate r, Y would be a constant. Further, the Eq. (3.26) shows that for 
p = 0, we also have M{r, t) = M(r) = constant. 

A. Collapse of a Physical Fluid 

We would emphasize that, for studying the collapse of a physical fluid, it is absolutely necessary to incorporate the 
radiation transport aspect in an organic fashion. The collapse equations were generalized to incorporate the presence 
of radiation by several authors |24|-|28|l . Following Misner [E4[ and Vaidya |2f| , we will first treat the radiation part of 
the stress energy tensor in the geometrical optics limit: 

E lk = qk l k k (3.28) 
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where q is both the energy density and the radiation flux in the proper frame and k % — (1; 1,0,0) is a null geodesic 
vector so that k % ki = 0. 

All one has to do now is to repeat the exercises for an adiabatic fluid outlined above by replacing the pure matter 
part of energy momentum tensor with the total one: 

T lk = (p + p)u l u k + pg lk + qk l k k (3.29) 

Then the new T® component of the field equation, upon integration, yields the new mass function: 

r-r r-r 

M(r,t)= AnR 2 dR(p + qv + q)= I dV[T{p + q) + qU] (3.30) 

Jo Jo 

Had we treated the radiation transport problem without assuming a simplified form of Eik and, on the other hand, 

in a most general manner, following Lindquist p8j, we would have obtained: 

M(r,t)= I dV[T{p + J)+HU] (3.31) 

Jo 

where 

J = E = E l UiUk = q = comoving energy density (3.32) 

and 

H = E° = average radial flux (3.33) 

This definition of M may be physically interpreted in the following way: while (p + q) is the locally measured energy 
density of matter and radiation, T(p + q) is the same sensed by S^ (r < 1). Here, the radiation part may be also 
explained in terms of " gravitational red -shift" . And the term HU may be interpreted as the Dopplcr shifted flux seen 
by Soo [gj,g5|,g8|. Although, the collapse equations, in general will change for such a general treatment of radiation 
transport, the generic constraint equation involving T incorporates this new definition of M and remains unchanged: 

r 2 = 1 + U 2 - ^- (3.34) 

R 

In the following, we list the other major collapse equations for the simplified form of E lk only: 

D T M = AttR 2 [T(p + q) + Uq] (3.35) 

D t M = -ATtR 2 pU - L(U + T) (3.36) 



where the comoving luminosity is 



L = AirR 2 q (3.39) 
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Even if there is no question of a strict exact solution (numerical or analytical) for such a fluid, it is believed by 
practically all the authors that a physical fluid will necessarily collapse to a singularity in a finite proper time; and the 
debate hinges on whether the singularity would be a BH or a naked one. The modern conviction in the inevitability of 
the occurrence of spacetime singularities in a general gravitational collapse of a sufficiently massive configuration stems 
on the strength of singularity theorems || . Probably, the first singularity theorem, in the context of spherical collapse, 
was presented by Penrose J29| where it was explicitly shown that once a trapped surface is formed, 2GM(r, t)/R > 1, 
the collapse to the central singularity is unavoidable. Since then many authors like Hawking, Geroch, Ellis, including 
Penrose himself, have proposed various forms of singularity theorems J3[ . In the next section, we shall show that the 
most innocuous assumption behind the singularity theorems, namely, the assumption that 

(1) the manifold should contain a trapped surface either in the past or future, 

is actually not obeyed by the collapse equations. 

It is clear that, if we are able to show that trapped surfaces are not formed even for the most idealized case of a 
nonrotating perfectly spherical perfect fluid not having any resistive agent like a strong magnetic field, certainly trapped 
surfaces would not form in more complicated situations. 

However, even before we present our deivation, it may be pointed out that we are aware of atleast one review 
article by Senovilla |29] which specifically describes the possibilty that the final state of a gravitational collapse may 
be singularity free (see subsection 7.2 of this reference). 

IV. PHYSICAL SPEED V 

It is of utmost importance to be able to properly define the quantity, V, which is the speed of the test particle or 
the fluid clement measured by a given static observer in a certain coordinate system. For a static field such as one 
considered here, the physical velocity in a generic coordinate system is given by using Eq.(88.10) or the the preceding 
unnumbered lowermost equation of pp. 250 of Landau and Lifshitz |[15| . In particular for the radial case, where, 
Va = Ve = 0, we have 



v _ dl _ V-9ndxi , 41 > 

dr y~g~oodx 

Here dxi/dxg is to be interpreted as the total or "convective" derivative, and is, in general, non-zero even if X\ is a 

comoving coordinate. In fact, except for Landau & Lifshitz, no other standard textbook on GTR seems to contain 

this discussion on the definition of V in detail, and, many experts on GTR also seem to be confused about this 

important aspect. One must note that, it is this V defined by Landaufe Lifshitz which appears in the Local 

Lorentz transformations. 

We have already discussed in Sec. Ill, the meaning of a dynamical derivative dr/dt in the COF. While in a 

Newtonian hydrodynamic problem, dr/dt defines the fluid speed, in the GTR case, following Landau & Lifshitz, the 

3-speed will be 



v _dl__ ^dr _ y-g rr dr 
dr Adt ^goodt 

It may be recalled again that while U(T) has the nature of a partial derivative with respect to t(r), it has the nature 
of a total derivative with respect to t(1). The reader can convince himself about this by going through references 
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p2| , [23|| . This is just similar to the fact while the fluid 4- velocity u l = ^j = ^- has a partial derivative nature with 
respect to t it is actually a total derivative of x l with respect to the proper time along the worldline of the fluid p8| . 
Then we note that the previously defined U and T are interlinked as 

U = TV (4.3) 

Here it may be also mentioned that the 3- velocity for the same fluid element may be different in different coordinate 
system. For instance, the 3- velocity measured in a Schwarzschild field would be 

eX/2dR ,A AS 

VscH - -^ (4.4) 

And, in general, Vsch 7^ T^(comoving). In particular, if the (vacuum external) Schwarzschild coordinates really, suffer 
from a coordinate singularity at R = R g , it is likely that Vs c h = 1 at R = R g and one would have Vs c h > 1 for R < R g 
(in the presence of matter-energy vacuum Schwarzschild singularity is irrelevant). Here one may argue that at the 
coordinate singularity, the Schwarzschild coordinates lose their meaning, and any result obtained by extending them 
further is devoid of physical meaning. However, since the comoving coordinates are by definition singularity free and 
one can always read off the time recorded by a comoving clock, we do not expect V > 1 in any truly non-singular 
domain of spacetime. 

V. THE CENTRAL PROOF: HEART OF THIS PAPER 

For purely radial motions, one may ignore the angular part of the metric to write: 

ds 2 = goodt 2 + g rr dr 2 (5-1) 

Again, by definition, the worldlines of photons or material particles are null or time like, i.e., ds 2 > 0, so that 

(grrdr 1 ' 



,9oo 



f 



\ goodt 2 



> (5.2) 



ffoo(l - V 2 ) > (5.3) 

Since in the singularity free COF, goo > 0, we find that j~ 2 = 1 — V 2 > 0, where 7 is the Lorentz factor of the fluid. 
This simple result expresses the fundamental fact that the 3-speed of photons or material particles cannot exceed the 
speed of light (as long as goo > 0). 

Further, we found in the Section III that there exists a quantity T 2 

and which is positive definite if so is A 2 = —g rr - 

Now, by substituting the above relationship (4.3) into the right hand side of another global constraint (3.22) or 
(3.34), we find that 
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Now by transposing, we obtain 



This equation may be rewritten as 



r 2 = i + r 2 y 2 - 2GM ^> (5.5) 

R{r, t) 



^-^- 2 -Wr <"> 



^ = 1- D ;;;; ( 5 - 7 ) 



r 2 _ 2GM(r,t) 

7 _1 ~ ^M) 

and this beautiful equation may be termed as the "master equation" for spherical gravitational evolution of a system 
of a fixed number of baryons. 

Since both T 2 and j 2 are positive definite, we find that the left side of Eq. (5.7) is positive definite; and so must 
be the right hand side of the same equation. Thus we obtain the most fundamental constraint for the GTR collapse 
(or expansion) problem, in an unbelieveably simple manner, as 

This shows that trapped surfaces do not form. 

Recall that it is believed that the energy of an isolated body cannot be negative [pl| . From physical view point, 
a negative value of Mb could imply repulsive gravity and hence is not acceptable. When we accept this theorem(s), 
we find that the fundamental constraint demands that if the collapse happens to proceed upto R — > 0, i.e., upto the 
central singularity, we must have 

M(r,i)->0; R^O (5.9) 

Remember here that the quantity M = mN (which is the baryonic mass of the star, if there arc no antibaryons) is 
conserved as Mf — > 0. Physically, the M = state may result when the negative gravitational energy exactly cancels 
the internal energy, the baryonic mass energy Mq and any other energy, and which is possible in the limit p — ► oo and 

p — ► oo. 

A. Singularity in Comoving Coordinates? 

Although, comoving coordinates, by definition, do not involve any singularity unlike external Schwarzschild coor- 
dinates, in a desperate attempt to ignore this foregoing small derivation, some readers might insist that there could 
be a coordinate singularity somewhere so that goo could be negative in Eq.(5.3). If so, (1 — V 2 ) would be negative 
implying that V > 1! And it might appear that, in such a case Eq. (5.6) would lead to 2GM JR > 1. First we want to 
emphasize that while the (External) vacuum Schwarzschild coordinate system might display such an anomaly, i.e, the 
3-speed of a free falling particle, measured in terms of i?and T would indeed appear to exceed the speed of light once 
it is inside the Event Horizon of a finite mass BH (if it would be allowed by GTR), the matter filled internal region 
would not display any coordinate singularity because the vacuum Sc. metric coefficients are completely irrelevant 
there. In fact the entire paper of Oppenheimer and Snyder B is aimed at finding the metric coefficients for the matter 
filled internal Schwarzschild coordinate system all the way upto R = 0, inside the supposed event horizon R < R g . In 
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particular, as emphasized again and again, the comoving coordinates are by definition singularity free. Even if one 
accepts this incorrect possibility for a moment, our eventual result survives such incorrect thinking in the following 
way. 

Since, the determinant of the metric, g = R A sin 9 goo g rr is always negative jTBj in all cases of coordinate singularity, 
when, goo < 0, we would have g rr > and A 2 < 0. Then it follows from Eq.(5.3) that T 2 < so that the L.H.S. of 
Eq.(5.6) would be again positive. Hence the R.H.S. of the same too must be positive. And we get back Eq.(5.8), and 
thus, it follows, in a most general fashion, that 

2GM ^ < 1 (5.10) 

R 

Hence, even if we, for the sake of arguments, accommodate the possibility (though quite incorrect) that the comoving 
coordinates might develop a coordinate singularity, our eventual physical result that trapped surfaces do not form 
remain unchanged. 

B. Previous Hints for M=0 Result 

While considering, the purely static GTR equilibrium configurations of dust, Harrison et al. |32] discussed long ago 
that spherical gravitational collapse should come to a decisive end with Mf = M* — 0, and, in fact, this understanding 
was formulated as a "Theorem" 

"THEOREM 23. Provided that matter does not undergo collapse at the microscopic level at any stage of com- 
pression, then, -regardless of all features of the equation of state - there exists for each fixed number of baryons A 
a "gravitationally collapsed configuration" , in which the mass-energy M* as sensed externally is zero." (Emp. by 
author) . 

In a somewhat more realistic way Zeldovich and Novikov Q| discussed the possibility of having an ultracompact 
configuration of degenerate fermions obeying the EOS p = e/3 with M — > and mentioned the possibility of having 
a machine for which Q —>■ MiC 2 . 

It is widely believed that Chandrasekhar's discovery that White Dwarfs (WD) can have a maximum mass set 
the stage for having a gravitational singular state with finite mass. The hydrostatic equilibrium of WDs can be 
approximately described by Newtonian poly tropes |34(| for which one has R oc pi , where p c is the central 

density of the polytrope having an index n. It shows that, for a singular state i.e., for p c — > oo, one must have 
R — ► for n > 1; and Chandrasekhar's limiting WD indeed has a zero radius [|34]]. On the other hand, the mass of 
the configuration M <x p c "• And unless n = 3, M — ► oo for the singular state. One obtains such a result for 

Newtonian polytropes because they are really not meant to handle real gravitational singularities. Fortunately, in the 
low density regime, when the baryons are nonrelativistic and only electrons are ultrarelativistic, the EOS is p — ► e/3 
and the corresponding n — > 3. Then one obtains a finite value of M c h - the Chandrasekhar mass. 

Now when we apply theory of polytropes for a case where the pressure is supplied by the baryons and not only by 
electrons, we must consider GTR polytropes of Tooper |3j|. It can be easily verified from Eq. (2.24) of this paper 
p5[ that in the limit p c — ► oo, the scale size of GTR polytropes A^ 1 — > 0. Further Eqs. (2.15) and (4.7) of the same 
paper tell that M ex K n > 2 ex p c — > for p c — ► oo. Thus, a proper GTR extension of Chandrasekhar's work would 
not lead to a BH of finite mass, but, on the other hand, to a singular state with M — > 0. 
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In a different context, it has been argued that naked singularities produced in spherical collapse must have Mf = 

s 

VI. PROBABLE REGIMES OF CONFUSION 

Although, the Positive Energy Theorems |3l| probe whether Mb cannot only be zero but even negative, and although 
many of the so-called naked singularity solutions correspond to zero gravitational mass J36[ , we can foresee that many 
readers would have difficulty in accepting our result. And, instinctively, there may be a tendency to reject this work 
on the basis of tangential and vague reasons. And though, we have taken great care in developing several ideas, in 
the face of the likely strong revulsion, several genuine or apparent confusions may creep up: 

A. Baryonic and Gravitational Mass 

For some readers it might appear that a M = state corresponds to zero baryon number N — 0. Very clearly, the 
reader, because of instinctive Newtonian notion, in such a case would incorrectly equate the gravitational mass with 
the baryonic mass: M = Mq = mN (incorrect). 

On the other hand, the gravitational mass of an isolated system is just the aggregate of all kinds of energy associated 
with it, and for any bound system, necessarily M < Mq. In particular, for the sake of illustration, we may recall that 
in the weak gravity regime, we would have 

M = Mq + E g + E m + E kineUc (6.1) 

Here the gravitational energy term is always negative (even if M < 0) and nonlinear. In the weak gravity regime 
it is ~ ~GM 2 /R, and as collapse proceeds, the grip of gravity becomes tighter, and this is effected by the non linear 
nature of E g . As a result, the value of M, in general, steadily decreases in any gravitational collapse, and, it is a 
natural consequence that if we have a continued collapse, the value of M will hurtle downward and the system would 
try to seek a state of "lowest energy" . In GTR, i.e., in Nature, the lowest energy corresponds to M = and not to its 
Newtonian counterpart En = E g + Ei n + EkineUc — (incorrect). Thus, if we remove the possibility of the occurrence 
of a repulsive gravity (negative M), then the bottom of the pit would be at Mf = 0. At this state, both E g and Ei n 
would be infinite but of opposite sign and separated by a finite gap Mq much like what happens in a renormalized 
Quantum Field Theory. 

There could be another confusion here as to how can | E g | be infinite when Mf = 0. This depends on how fast the 
value of M — > with respect to R — » and is perfectly allowed for a singular state. 

Nevertheless, it must be borne in mind that for an actual strong gravity case all such contributions shown in the 
foregoing equation intermingle with one another in a non-linear and inseperable manner. 

B. Principle of Equivalence (POE) 

Even though there are many published results suggesting M = in connection with naked singularities, our work 
might be singled out with the plea that a M = result violates POE. We repeat once again that, POE only says 
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that the local nongravitational laws of physics are the same as the corresponding laws in STR. For example, this 
would mean that the Stefan- Boltzman law which tells that the emissivity of a black body surface is oc T 4 , remains 
unchanged. POE does not impose any limit on the value of T itself and hence on the total amount of radiation emitted 
from the black body surface. POE does not say that only a certain percentage of the initial total mass energy Mi can 
be radiated in the process, POE has got nothing to do with either the imposition of any additional local constraint 
(such as a maximum value of T) or any global issues. 

If one would invoke POE to debar phenomenon which are not understandable in Newtonian notions (like M = Mq, 
incorrectly) GTR itself is to be discarded. With such a viewpoint, all work on Positive Energy Theorems are to be 
considered as redundant and unnecessary because in STR, the mass-energy of a system which was positive to start 
with can never be negative. 

C. Matter - Antimatter Annihilation ? 

In STR, there is no gravity and hence there is no Kelvin- Hclhmoltz process, neither could there be any real finite 
material body held together by any long range force ( a plasma has to be confined by external electromagnetic fields). 
And there could be a naive idea that the entire initial mass energy may be radiated only if there are processes like 
e + e~ — ► 27. If this is envisaged as the only way to generate radiation (in this case photons), it must be remembered 
that such a thing refers to systems having total lepton number or total baryon number as zero. For matter consisting 
of a definite baryon number and lepton number there cannot be any energy extraction by this process. Yet such 
matter radiates because of normal electromagnetic processes like Bremsstralung or Compton processes, or by nuclear 
processes like p + p — > it — > 27. Actually at very high densities and temperatures, in astrophysical scenarios, energy 
is liberated by the so-called URCA or weak interaction processes involving emission of vv. Whatever be the process, 
if the global Kelvin- Helmholtz process heats up the matter to sufficiently high temperature near the singularity (to 
which everybody agrees), the center of mass energy of the colliding particles, like, electrons, protons, neutrons, quarks 
or whatever it may be, will be accordingly high enough. And in this limit, for an individual collision, the colliding 
particles can radiate not only an energy equal to their rest mass but any amount higher than this. The easiest example 
would be that an e~ — e + collider can generate particles (photons, neutrinos, quarks etc.) much heavier than 0.5MeV. 
And it should be also remembered that when we say that the entire MiC 2 may be radiated, we do not mean that this 
happens in a flash as is the case for matter-antimatter annihilation. On the other hand, in gravitational collapse, it 
is the integrated radiation over the entire history of the process we are concerned with. 

VII. DUST COLLAPSE 

There is no way we can ever think of exactly solving even the adiabatic collpase equations for a real fluid, i.e., one 
having pressure. Further, this idea of adiabaticity would break down as soon as the fluid starts to contract because 
of the Kelvin- rlelmholtz energy liberation. Even if we consider the fluid to be degenerate and at T = to start with, 
gravitational contraction would keep on heating it up unless it acquires an effective adiabatic index j t — 4/3. On the 
other hand, we may feign to ignore the role of any temperature in the fluid by artificially assuming a polytropic EOS, 
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p ex p 7t even when the gas is non- degenerate. But the value of 74 will keep on evolving and it is not possible to find 
any unique solution for the entire range of p and p even by any numerical means. Depending on the inevitable hidden 
assumptions made, it may be possible to obtain various solutions and none of which may have to do much with the 
actual complicated physics of atomic and nuclear matter at arbitrary high density and pressure (z s — 00). And, the 
only way, one may hope to obtain an exact or near exact solution, at the cost of the actual thermodynamics, is to do 
away with the EOS, i,e., to set p = even when p — > 00!! Even then, there could be exact analytical solution only 
when the dust ball has uniform density. And this problem was first (apparently) solved by Oppenheimer and Snyder 
(OS) ||]. It may be mentioned that recently the problem of the occurrence of the final state in a GTR collapse, 
with the inclusion of pressure gradients, have been considered by Cooperstock et al. 0. These authors have correctly 
pointed out that at very high densities the effect of pressure gradient must be included and for the likely occurrence 
of negative pressure gradients, there may may not be BH formation at all. 

OS initially worked in the COF, but, then, to match the internal solutions with the external ones, eventually 
shifted to the (non comoving) Sch. frame involving R, T. Without giving the details of the actual mathematical 
manipulations, we shall simply present their key equations. By matching the internal solutions with the exterior ones 
they obtained a general form of the metric coefficients and also a relation between T and R which is valid for the 
entire range of 00 > R > : 

grT^e^iidT/drfil-U 2 )}- 1 (7.1) 

-9RR = e x = {l-U 2 )- 1 (7.2) 



and, 



where 



T = - 3 nJ"(rr _ fiVy' 2 ) - 2JW" + «.» i" C^ P-3) 



»- 5 [<"•»>' -^^ < 7 - 4 ' 

It is the above Eq. (7.3) which corresponds to Eq. (36) in the OS paper. OS also showed that the relation between 
T and r is determined by 

Ft + r 3/2 = R 3/2 (7.5) 

where, 

F=-(3/2)RH 2 (r/r b ) 2 ; r < r b (7.6) 

So, for the outer boundary, we have 



2 r 3/ 2 _ _R3/2 

T = — - 



3 [r/nYR]' 2 
According to OS, in the limit of large T, one can write 



(7.7) 
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The last term of the above equation contains a typographical error, and, in general, this equation missed a numerical 
factor of 4. The corrected form should be 



T~-R gb In 1 



2 
^6 



" ^-^n (r.., 



4i? s6 V 2r fc 3/2 



From the foregoing equation, they concluded that, "for a fixed value of r as T tends toward infinity, r tends to a finite 
limit, which increases with r". 

— 1/2 

It follows that the supposed finite limit for r ex R , . OS then found the explicit expressions for the internal Sch. 
metric coefficients in the limit of T — > oo. However, they missed a numerical factor of 4. And though this numerical 
factor does not change the qualitative behavior of the solutions, in the following, we present the corrected relevant 
expressions 

er X = 1 - (r/r b f Ue- T / R ^ + ± [3 - (r/r b f] } (7.10) 



and 



X-2T/R gb ( 1 

e -T/* sb + 1 [3 _ (r/rb) 2| . (? n) 



4 [ 2 

Note that these equations were obtained by eliminating R, and the T = oo limit covers the entire R < i? g f, range. 
Also note that the comoving coordinates r and rt are fixed, and when there is a complete collapse to a physical point 
at R = 0, e A ?tjms£ Mow) «p irrespective of the value of r. But this does not happen for the OS solution for an 
internal point r < r b \ 

For instance if R g b > 0, for rjr^ = 0.5, we obtain e A = 15/13 at T = oo, when we should have obtained e A = oo. 
And for the central point r = (which also corresponds to R = 0), we find that after the collapse (T = oo), we have 
e A = 1 when we should have had e A = oo! And OS, somewhat casually, noted this: "For A tends to a finite limit for 
R < Rgb as T approaches infinity, and for Rb = R g b tends to infinity. Also for R < R g b, v tends to minus infinity." 

However, unfortunately, OS did not bother to ponder on the genesis of this completely unphysical aspect of their 
solution. This, on the other hand, is a definite signature that there is a severe problem in the foundation of 
this problem. 

A. True Solution of the O-S Problem 

We note that OS completely overlooked the most important feature of Eq. (7.3) (their Eq. 36), that in view of the 
presence of the T ~ In v 1 / 2 +1 l term, in order that T is definable at all, one must have y > 1 

For an insight into the problem, we first focus attention on the outermost layer where yb = Rb/R g b, so that the 
above condition becomes 

2GM b /R b c 2 < 1 (7.12) 
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Thus a careful analysis of the GTR homogeneous dust problem as enunciated by OS themselves actually tell that 
trapped surfaces cannot be formed even though one is free to chase the limit R — ► R gb . This means that, the final 
gravitational mass of the configuration is 

M f (R = 0)=M g (R b = R gb )=0 (7.13) 

But then, for a dust or any adiabatically evolving fluid Mi = Mf — constant. Therefore, we must have Mi — too. 
And for a finite value of R, this is possible only if p = 0. But for a dust p = po = mn and therefore, we have n = 0. 
Finally, the total number of the baryons in the configuration N = 0. From, a purely mathematical view point, the 
N = limit can be described as 

r = r b -* 0; r/r b -» 1 (7.14) 

Further note that if we really assume R gb 7^ 0, the second term r b R/R gb r of Eq.(7.4) can be made arbitrarily 
small as the collapse proceeds to R — > 0. Remember here that r and r b are comoving coordinates and are fixed by 
definition. So for any interior region r separated from the boundary by a finite amount r < r bl y becomes negative 
in contravention of Eq.(7.3) if R gb 7^ 0: 

y -+ ^ [( r /n) 2 - 1] < as R -► (7.15) 

This is alleviated if either or both of the two following conditions are satisfied : (i) r = r b , as derived above or (ii) 
Rgb = 0, which again leads to the previous condition. 

Thus had O-S carefully noted this simple point, they would probably not have proceeded with the rest part of their 
paper which hints at the formation of a finite mass BH in a completely erroneous manner. And mathematically, the 
Eqs. (7.10) and (7.11), in a self-consistent manner degenerate to definite limiting forms: 

e- x w 1 - Ue- T/R " h + 1\ -» 0; e A -* 00, (7.16) 

e v a e — A { e - T / R s> + l} - (7.17) 



B. Common Perception About Formation of Horizon 

There is a widespread idea that atleast for a dust collapse, irrespective of our above explicit proof, the formation 
of an Event Horizon is inevitable. One assumes here a shell of dust particles to be either at rest or in equilibrium 
at t = 0. Then suddenly the dust is envisaged to lose its state of rest and is allowed to collapse. And since a dust 
does not radiate its gravitational mass remains fixed and it is expected to reach its horizon R g — 2GM/c 2 in a finite 
proper time on the basis of the OS solution (which ignored the y > 1 condition). 

For a resolution of this puzzle, first one has to appreciate that if a dust ball is ever at rest without the support of 
external mass energy, its mass must be zero. To see this simply consider the Oppenheimer-Volkoff equation fl37fl for 
hydrostatic balance: 

dp p + p(0) 



dR R(R - 2GM) 
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(AirpR 6 + 2GM) (7.18) 



If p = 0, the above equation yields p — too |3£| (Private Comm., A.K. Raychowdhury) . From the view point of 
thermodynamics, any physically meaningful EOS will yield p — if p = irrespective of whether the fluid is in 
hydrostatic equilibrium or not. And when p = everywhere, it can be shown that the proper time for collapse to 
both the horizon and the central singularity is oo. In particular, the latter is given by 

where p(0) is the density of the dust when it was at rest; U(0) = 0. 

For a finite value of i?&(0) = Rq, one has M = for p = 0. And if Ro — oo to make M > 0, one would again 
find the proper time to reach the horizon to be oo. Now the reader is requested to carefully go through the following 
subtle points for a resolution of this puzzle. 

• One can imagine the dust ball to be in a state of rest by either unconsciously assuming the presence of (1) some 
finite pressure gradient forces or (2) some external source of mass energy, like, friction forces or biological forces at 
t = 0. The latter possibilty eventually boils down to the presence of some external electromagnetic force field in the 
problem at t = 0. 

• Thus at t = we do not have a dust ball, but, on the other hand, we have either (1) a physical fluid or (2) a 
system of particles interacting with some external (electromagnetic) sources of mass-energy 

• If the fluid loses its hydrostatic balance, surely, it would start collapsing at t = 0, but it does not mean that 
the collapse would be a free fall. In other words the relevant initial condition of the problem that there is a "sudden 
collapse" can be realized, but, what cannot be realized is the idea that p = finite for t < but again p = for t > 0. 
This assumption that pressure vanishes instantly violates causality and is not allowed by GTR. 

• In the latter case too, the external electromagnetic field must decay over a finite time scale, howsoever, small. 
And when this finite time scale removal of the external electromagnetic field will be implemented, the problem would 
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be different from the ideal "dust collpse" problem. The OS result r oc M i will cease to be valid in such a case. 
And we must not let loose our intuitive Newtonian concepts to determine the fate of such a problem by overrdiding the 
general constraint 2GM/R < 1. In particular, here R must not be confused as the proper length along the worldline. 

• A pure dust is however allowed to collapse from a state of rest by starting from oo. But once it statrs from oo, it 
would never reach the would be EH supposed to be located at a finite radial distance in a finite comoving time. This 
would be in agreement with our contention. 

VIII. KRUSKAL COORDINATES 

Since we have already shown in a most general manner that because of an inherent Global Constraint, the Einstein 
equations dictate that if the collapse proceeds upto R = 0, the gravitational mass of the singularity would be M = 0, 
it is clear that, if there would be any BH formed by gravitational collapse, its mass would be zero. To show the 
consistency of our result, in the following, we assume first the existence of a finite mass Sch. BH, and then, try to 
verify whether a finite value of M > is allowed or not. Both the exterior and interior spacetime of a BH is believed 



to be described by the Kruskal coordinates |39|. For the exterior region, we have (Sector I and III): 



u = A(.R)cosh— -; v = MR)swh- T7 ) R>2M (8.1) 

AM AM ~ 
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where 

h{R) = ±(^-lj' 2 e R/AM (8.2) 

where the +ve sign refers to "our universe" and the -ve sign refers to the "other universe" implied by Kruskal diagram. 
It would be profitable to note that 

d/i ±R ( R \ _1 2 r/4m /oo. 

dR ~ m* \2M ~ 1 ) e (8 - 3) 

And for the region interior to the horizon (Sector II and IV), we have 

u = f 2 (R)smh^; v = f 2 (R) cosh ^ ; R < 2M (8.4) 

where 

MR) = ± ( 1 -t^) 12 z R/4M = V=ih(R) (8.5) 

where, again, the +vc sign refers to "our universe" and the -ve sign refers to the "other universe". It is found that, 

-1/2 



#2 = t-R L R\ 



dR 8M 2 V 2MJ [ ! 

Given our adopted signature of spacctimc (—2), in terms of u and v, the metric for the entire spacetime is 

ds 2 = : ^Ml e - R / 2M (dv 2 - du 2 ) - R 2 (d9 2 + dej? sin 2 9) (8.7) 

R 

The metric coefficients are regular everywhere except at the intrinsic singularity R = 0, as is expected. Note that, the 
angular part of the metric remains unchanged by such transformations and R(u, v) continues to signal its intrinsic 
spacelike nature. In either region we have 

u 2 -v 2 = (— - l) e R ' 2M (8.8) 

so that 

u 2 -v 2 ^0; u 2 ^v 2 : R = 2M (8.9) 

In the External region, we also have 

u 

- = coth(T/4M) (8.10) 

v 

For studying the Kruskal dynamics, it would be useful to briefly recall the dynamics of a test particle in the Sch. 
coordinates because the Kruskal coordinates arc built by using R and T . 



A. RADIAL GEODESIC IN SCHWARZSCHILD COORDINATE 

For a radial geodesic, we have 

ds 2 = dT 2 (l - 2M/R) - (1 - 2M/R)- 1 dR 2 (8.11) 
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In terms of the conserved energy per unit rest mass 



it is possible to find that @@ 



jrj~\ 

E = —(1-2M/R) (8.12) 

UT 



d ^ = - l —fl^[E*-(l-2M,R)]^ (8.13) 



If E — > oo, as it happens for a photon, for arbitrary value of 2M/R, the foregoing equation attains a form 

g->-(l-2M/fl) (8.14) 

Interestingly, even when E is finite, Eq. (8.13) attains the same form if R — ► 2M. 



B. Radial Kruskal Geodesic 

We would like to explicitly verify whether the (radial) geodesies of material particles are indeed timelike at the EH 
which they must be if this idea of a finite mass Schwarzschild BH is physically correct. First we focus attention on 
the region R > 2M and differentiate Eq.(8.1) to see 

du du dudT df , T f . T dT . r . 

lR = M + dfdR = dR C ° Sh lM + lM Smh lMdR (8 ' 15) 

Now by using Eqs. (8.1-3) in the above equation, we find that 



du ru , „ , „ i v dT 

dR = 8M^ 2M -^ + 4MdR 



^ = ^W2M-l)- 1 + ^7^; R>2M (8.16) 



Similarly, we also find that 



dv rv ,„,„,, JN _i u dT 



— = ^(R/2M -I)- 1 + ; R>2M (8.17) 

dR 8M 2V ' ' U'IdR' ~ v ' 

By using Eq. (8.13) in the two foregoing equations and then by dividing (8.16) by (8.17), we obtain 

vE 



du y/E 2 -l+2M/R 

dv ~ ~ ^ 



(8.18) 



^/E 2 -1+2AI/R 

Interestingly, we have verified that one would obtain this same equation for du/ dv even for Sectors II and IV (R < 2M). 
Since u and v are expected to be differentiable smooth continuous functions everywhere except at R = 0, and also 
since the "other universe" is a mirror image of "our universe" , we expect that the value of du/dv for any given R 
must be the same, except for a probable difference in the signature, in both the universes. For the "our universe" , it 
is found that, we have 

UH = v H (8.19) 

where "H" refers to the value on the event horizon R = 2M . And it can be found that if one approaches the horizon 
from Sectors II or IV (other universe), one would have 

u H = -vh (8.20) 
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In this case, we find 

du uh — vh 2uh 



= -1; R = 2M (8.21) 



dv vh — uh 2vh 
so that, we have du 2 = dv 2 at R = 2M. Then, we promptly find that for a radial geodesic 

ds 2 = 16M 2 e- 1 (du 2 -du 2 ) = 0; R = 2M (8.22) 

This implies that although the metric coefficients can be made to appear regular, the radial geodesic of a material 
particle becomes null at the event horizon of a finite mass BH in contravention of the basic premises of GTR! 

And since, now, we cannot blame the coordinate system to be faulty for this occurrence, the only way we can 
explain this result is that the Event Horizon itself corresponds to the physical singularity or, in other words, the mass 
of the Schwarzschild BHS M = 0. 

What would happen if we would approach the EH from sector I to III? In this case, by taking appropriate limit, it 
can be shown that 

du = (e/ug) - (u H /E 2 ) 

dv -(e/u H ) + (uh/E 2 ) [ ■ ' 

Here the magnitude of du/dv appears to be a function of E if uh = Vh is finite. In particular, for E = I, again 
du/dv = — 1 irrespective of the value of ujj and again ds 2 =0on the EH. Note that, as long as ds 2 > (time like), its 
value may be allowed to be a function of the initial conditions, but once its value is (null) for any initial condition, 
for the sake of consistency, its value must be for any other initial condition too. This is similar to the following 
situation: As long as the value of V < 1, its value may be different for different observers and initial conditions. But 
once any observer finds V — ► 1, all other local Lorentz observers would also find V — ► 1. Coming back to the case of 
du/dv, in order that du/dv = — 1 at the EH for all the observers, we must have uh — oo, which demands M = 0. 

And then, the entire conundrum of "Schwarzschild singularity" , "swapping of spatial and temporal characters by 
R and T inside the event horizon {when the angular part of all metrics suggest that R has a spacelike character even 
within the horizon), "White Holes" and "Other Universes" associated with the full Kruskal diagram get resolved. 

C. Physical speed at the Horizon 

As mentioned earlier, the physical velocity of a particle under the influence of a static gravitation field, as measured 
by a certain local observer can be found by using Eq.(88.10) or the the preceding unnumbered lowermost equation of 
pp. 250 of Landau and Lifshitz JLq]: 



V = - = ^sT ldxi (8.24) 

dr y/g^dxo 

Thus for the Kruskal case, the radial speed of free fall is 

Vk = P (8.25) 

dv 

While extending this idea of a "locally measured" 3-speed one point is to be borne in mind. As emphasized by Landau 
and Lifshitz [|15( , there cannot be any "static observer" at the Event Horizon. However, we can conceive that there 
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is a static observer at R = 2M + e, where e — > 0. In other words, we may only study the limiting behaviour of V at 
the EH. Thus, our result shows that the free fall speed at the EH, V — > 1, and this is not allowed by GTR unless 
R g b = 2M = 0. Now we explain why Vk — ► 1 at the EH for any coordinate system, Kruskal or Lemaitre or anything 
else. Let the speed of the static other observer be Vsch-o with respect to the Schwarzschild observer. By principle of 
equivalence, we can invoke special theory of relativity locally. Then the free fall speed of the material particle with 
respect to the other static observer will be 

T/ Vsch ± Vsch-O to oc x 

v - TTv — v ^ b > 

1 ± v Sch VSch-O 



We also find that for the Sch. metric |(T^,|T^] , we have 



Vsch — 



1 / 2AA 



.27) 



Note that for a photon, E = oo, and the above Eq. correctly yields V v hoton = 1 anywhere. And for a material particle 
having E > 0, it shows that, again Vsch — ► 1 as R — > 2M. Consequently, at the Event Horizon, Eq.(8.25) would 
always yield | V |— > 1 too. The value ofV can change in various coordinates only as long as it is subluminous to all 
observers. 

D. Back to Schwarzschild Coordinate System 

We found that ds 2 (u,v) — at R — 2M . But actually ds 2 is invariant under coordinate transformation, and, there 
it should be obvious that ds 2 (R,T) — too at R — 2M. And it is easy to verify that it is indeed so. We found in 
Eq.(8.14) that either for a photon anywhere or for a material particle at R — > 2M, 

{dR/dT) 2 -»• (1 - 2M/R) 2 (8.28) 

Infact this above condition is called "null geodesic" one, and we find that, when this is satisfied, we have 

dR 2 

1 _ 2M/R - C 1 " 2A//i?)dT 2 = dz 2 (8.29) 

By using this Eq. in Eq.(8.11), we find that for radial geodesies, either for a photon anywhere or for a material 
particle at R = 2M, we have 

ds 2 (R, T) -> dz 2 -dz 2 = (8.30) 

Normally one would ignore this result as a reflection of the" coordinate singularity". But, even then, the value of 
ds 2 is an invariant, and its value must be same in all coordinates, as we have already verified for the Event Horizon. 
Sadly, like a proverbial "darkness beneath the lamp" this simple point was not seriously taken note of in the past! 

IX. SUMMARY AND CONCLUSIONS 

We found that for the continued collapse of any perfect fluid possessing arbitrary EOS and radiation transport 
properties, a proper amalgamation of the inherent global constraints arising because of the dependence of spatial 

25 



curvature like parameter (r) on the global mass-energy content, M, directly shows that no trapped surface is allowed 
by GTR. This result becomes independent of the details of the radiation transport properties because the integration 
of the (0,0)- component of the Einstein equation, yields the definition of M by absorbing all quantities like p, q, and H , 
in whichever fashion they may be present. Then it follows that if there is a continued collapse, the final gravitational 
mass of the configuration necessarily becomes zero. This Mf = state must not be confused as a vacuum state, on 
the other hand, the baryons and leptons are crushed to the singularity with an infinite negative gravitational energy 
Eg — > — oo. On the other hand, the positive internal energy is also infinite Ei n — > oo. However, in the present paper, 
we did not investigate whether this state corresponds to 2GMf/R < 1 or 2GMf / R = 1. In another work [J40[ , we find 
that, it is the latter limit which should be appropriate, i.e., the system keeps on radiating and tends to attain the 
state of a zero mass BH characterized by zero energy and entropy, the ultimate ground state of classical physics. Since 
the concept of a BH is intrinically meaningful only for M > 0, by borrowing a terminology from Cooperstock et al. 
|| , we may call this final state to be marginally naked. In fact the Eq. (17) of this paper |]8J considered the possibility 
that as R — ► 0, one may have M — > 2M — > 0. And it is precisely this result which we have obtained without making 
any kind of assumption. 

In this paper, we did not try to find the proper time required to attain this absolute classical singular ground state 
though we found that for the fictitious dust solutions r = oo. This question has, however, been explored elsewhere 
p0[ to find that, for a real fluid too, r = oo. This means that there is no incompleteness in the radial worldlines of the 
collapsing fluid particles inspite of R having a finite range (rb). Such a Non-Newtonian behavior is understandable in 
GTR because, it was found jfCJ] that although M keeps on decreasing, the curvature components ~ GM/R 3 ~ R~ 2 
and — g rr tend to blow up. As a result the 3-space gets stretched and stretched by the strong grip of gravity, or in 
other words, the proper distances eventually tend to blow up too. 

We also found that this inherent global constraint is also imprinted in the important work of Oppenheimer and 
Snyder g because in order that, at the boundary of the star, 

T~ln4^=ln^C±i (9.1) 

yV2_i [R b /R gb )V 2 -l v ; 

remains definable, one must have 2GMb/ Rb < 1. And then the central singulaity Rb — > could be reached only if 
Mb = 0, if the horizon coincides with the central singularity. Accordingly, the value of r g b = oo along with T = oo. 
In fact this result follows in a trivial fashion from Eq. (32) of their paper (our Eq.[8.10]) 

y^\[{r/nf-i] + r ^ (9.2) 

2 r Rgb 

where the parameter y which must be positive. But if R g b ^ 0, as R — ► 0, it is trivial to see that y actually 
becomes negative for r < rb- This shows that actually the horizon or any trapped surface in never allowed by the 
OS solution. And we noted that since OS did not incoroprate this intrinsic constraint in their eventual approximate 
expression for e A , it failed to blow up for the internal regions of the dust ball even when the collapse is 
complete. 

Despite having proved our results in a general fashion, we reconsidered the case of a finite mass Sch. BH described 
by Kruskal coordinates. We found that the radial geodesic of a material particle, which must be timelike at R — 2M, if 
indeed M > 0, actually becomes null. And this independently points out that M — 0. This simple fact independently 
asserts that there is no Event Horizon, no Schwarzschild Singularity, no T-rcgion, and the only singularity that might 
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have been present is the central singularity, and whose mass must be zero. Technically, one might view this central 
singularity as the Schwarzschild Singularity associated with a zero mass BH. Even then the existence of such a zero 
mass BH could be realized only if the collapse process could be complete in a finite proper time; but it actually takes 
infinite time : Nature abhors not only naked singularities but all singularities; and we find that only GTR may be 
having the mechanism of removing such singularities even at a classical level. This precise possibilty has recently 
been considered by Senovilla p9J] . And this happens because of the marriage between the physics and (spacetime) 
geometry. If somehow, one would try to build up a super concentrated energy density near a "point" , the space 
would get dynamically stretched by the gravity associated with the concentrated energy density and a singularity is 
avoided. But we emphasize that the present discussion does not rule out the likely occurrence of a singularity in the 
cosmological context. 

Consequently, several associated theoretical confusions like (i) whether the physically defined circumference coor- 
dinate R can, suddenly become a time-like coordinate, (ii) whether there could be White Holes freely spewing out 
matter and energy in the observable universe, and (iii) whether information can really be lost from the observable 
universe in violation of the quantum mechanics, which have plagued GTR in the present century, would be resolved, 
if the present work is correct. 

Finally, we appreciate the physical intuition of Einstein [tt2| and Landau p3] in not being able to accept the reality 
of Schrawzschild Singularity or any singularity in GTR. We also recall that Rosen p4| , in an unambiguous manner 
noted the impossible and unphysical nature of the T-region. 

"so that in this region R is timelike and T is spacelike. However, this is an impossible situation, for we have seen 
that R is defined in terms of the circumference of a circle so that R is spacelike, and we are therefore faced with 
a contradiction. We must conclude that the portion of space corresponding to R < 2M is non-physical. This is a 
situation which a coordinate transformation even one which removes a singularity can not change. What it means is 
that the surface R — 2M represents the boundary of physical space and should be regarded as an impenetrable barrier 
for particles and light rays." (Emp. by author). 

We have, in this paper, attempted to resolve all such paradoxes by showing that not only the R < 2M region 
unphysical, it does not exist or is not ever created. 

We have pointed out several instances when it was hinted or suggested that the GTR singularity may correspond to 
AI = state. Additionally we point out that the numerical studies of collapse of scalar fields suggest that it is possible 
to have BHs of M — |4j| . More importantly, the supersymmetric string theories find the existence of extremal BHs 
with charge Q — M, which for the chargeless case yields M = ]45j| . However, ironically, string theorists, at this 
moment, are guided by the erroneous notion that GTR yields BH with M > 0, and are struggling to wiggle out of 
this result by modifying the definition of event horizon into the socalled "stretched horizon" . 

Although, it might appear that astrophysics would be poorer in the absence of the mystique of BH, actually, it 
may be possible to envisage new varieties of stable or quasi- stable ultracompact compact objects of stellar mass or 
dynamically contracting super massive stars responsible for new gamut of astrophysical phenomenon. But there is 
one constraint imposed by GTR on such probable static ultracompact objects : if the compact object is assumed to 
be cold and in hydrostatic equilibrium, the surface redshift z s = T _1 — 1 < 2 lira, and it does not mean that there can 
not be any compact object beyond this limit, i.e., z s > 2. It only means that such high z s > 2 objects must be "hot" 
and dynamically contracting (remember the time to collapse to a singularity is oo). Also, in case there is a positive 
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(repulsive) cosmological constant, it may be possible to have more massive static compact objects. 

Also recall here that if rotation is taken into considerations, the value of Mov could be significantly raised. But 
to fully appreciate the question of likely existence of stellar mass BH candidates of masses as high as ~ WM@, we 
must keep an open mind with regard to our present day understanding of QCD. Even with reference to present state 
of knowledge of QCD, there could be compact objects with exotic EOS, where the masses could be ~ 10M Q or even 
higher |46| . These stars are called Q-stars (not the usual quark stars), and they could be much more compact than 
a canonical NS; for instance, a stable non rotating Q-star of mass 12Mq might have a radius of ~ 52 Km. This may 
be compared with the value of R g b ~ 36 Km of a supposed BH of same mass. 

In general, it is believed that, at sufficient high temperature, quark confinement may melt away. And the energy 
gained from the pairing of quarks and antiquarks of all colors which drive the chiral symmetry breaking may be 
overcome by the entropic advantage in letting the particles be free. At a very high T, therefore, asymptotically, free 
quarks, antiquarks and gluons should be liberated |47(| and provide new sources of pressure. There is already some 
evidences that at a temperature of ~ 150 MeV, there is a phase transition in hot nuclear matter and new degrees 
of freedom are suddenly liberated. It is such processes which may allow ultracompact objects to be in a stable or 
dynamic quasi-stable state. 

However, the above argument that there there may be static ultracompact objects having masses larger than 
canonical neutron stars does not at all mean that we are advocating that all the central compact objects hinted in 
X-ray binaries or Active Galaxies must be static. On the other hand, our work shows that there could be compact 
objects which would take infinite proper time to collapse a singular state because although, in the coordinate space 
the objects may be sinking, in the inner physical space (proper radius), they may be expanding! In other words, the 
final state suggested by our work (for a physical fluid) is the one with T — > 0, U — > 0, but with V — > c. Recall, when 
we say "radius" of a body we imply R, and thus the rate of observed radial contraction in such cases would be given 
by dR/dT — > 0. Yet, in internal space, the body may be contracting with a speed V — > c ! But, at a finite proper 
and lab time, when we detect a massive condensation or ECO, it is likely to have a value of R almost equal to its 
Schwarzschild radius and the value of goo on the surface could be exremely low but finite. For example it may have 
a value of 500 = 10~ 10 , 10 -11 or even less (but finite). Since, the corresponding lab speed of contraction dR/dT w 0, 
the ECO would look like a static and frozen object over time scale of years. Yet, its actual local speed of collapse 
(in the space of proper radial distance) V ~ (<7oo) — 1 dR/dT « c or also V <C c. We understand that it is extremely 
difficult to comprehend this because, intuitively, we tend to equate R as the measure of proper radial distance. 

When advection dominated accretion flow will interact with such an Eternally Collapsing Object (ECO), in case, 
indeed, 7«c, the collision process would emit insignificant radiation. However, if one assumes that the central object 
is static, one would expect large luminosities from the surface of the object. And the absence of such luminosities (in 
case V ~ c) would be interpreted as the "evidence for an event horizon" ! An ECO is struggling to attain the R = 0, 
<7oo = an d g r r — 00 state, but it would succeed to reach there only after the elapse of infinite proper time by which 
time its gravitational mass M — ► 0. 

Let us briefly recall the case of the recently discovered Cosmological Gamma Bursts like GRB990123 J48| and 
GRB980329 J49] which have been estimated to have radiated equivalent of Q 7 ~ 2M Q c 2 in the electromagnetic band 
alone under assumption of isotropy. The total energy radiated including the neutrino emission is expected to be 
atleast twice this amount Q ~ 4Af Q c 2 . And this is in agreement with the predicted and estimated true energy budget 
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of poweful Gamma Ray Bursts [Q . It is possible that many GRB afterglows are highly beamed and the actual energy 
budget could be considerably lower than estimated assuming isotropy. Note, while a BH cannot have any intrinsic 
magnetic field, an UCO/ECO could be highly magnetized and thus the latter is much more capable to explain likely 
beamed emission. In fact, the observed spectral break and rapid fading of the optical afterglow [pTJ in GRB990123 
has been interpreted in terms of beaming. But this interpretation is not fully satisfying because (i) the spectral break 
may be explained by spherical models fl52] and (ii) no linear polarization has been detected in the optical afterglow 
|J53| . In general, the long term afterglow observations of GRB 970508 and 971214, show that they are inconsistent 
with jet models Q. Further the afterglow of GRB 970402, 970616 and 98042588 too are consistent with isotropic 
models pa ]. Such energy emission can be explained as collapse of massive stellar cores, and is hardly possible if 
trapped surfaces really formed at values of Mf « Mj. On the other hand such phenomenon might be signalling the 
formation of new relativistic ultracompact objects or ECOs. 
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